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Abstract 

A simple model to explain numerically observed behaviour of chaotically varying stripes and 
square patterns in zero Prandtl number convection in Boussinesq fluid is presented. The nonlinear 
interaction of mutually perpendicular sets of wavy rolls, via higher order modes, may lead to a 
competition between the two sets of rolls. The appearance of square patterns is due to secondary 
forward Hopf bifurcation of a set of wavy rolls. 
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The study of low Prandtl number thermal convection [1-13] has been long motivated by 
its importance in geophysical and astrophysical problems. The hydrodynamical equations of 
thermal convection in Boussinesq fluids involve two types of nonlinearity. The first describes 
self-interaction v.Vv of the velocity field v, and the second nonlinearity v.V# results from 
the advection of the temperature fluctuation 9 by the velocity field. The nonlinearity v.V# 
may be neglected in the asymptotic limit of zero Prandtl number |rj. The linearly growing 
two-dimensional (2D) rolls then become exact solutions of nonlinear equations, if stress-free 
boundary conditions are considered. This makes this limit interesting even from purely the- 



oretical point. Thual |12j recently showed by a most general 3D direct numerical simulations 
(DNS) of zero P asymptotic equations that the solutions do not blow up. The comparison 
with full Boussinesq equations with both nonlinearities also reproduced zero P results. This 
DNS showed many interesting patterns including the possibility of square patterns in zero P 
Boussinesq fluids with stress-free boundaries for 1.05 < r{= R/R c ) < 1.7, where the critical 
Rayleigh number R c = 277r 4 /4. However, the mechanism of generation of square patterns 
in zero P convection remains unexplained. A nonlinear interaction between 2D rolls can- 



not generate either square or hexagonal patterns [ 1 1 for zero P asymptotic equations as 
growing straight 2D rolls are their exact solutions. The streamlines in DNS |]12| support 
this view. On the other hand, the mechanism of prevention of continuous growth of 2D 
rolls approximately 0.5% above the onset of convection was captured in a simple dynamical 
system ||13|| , which agreed well with the results of DNS in its validity range. This model 
suggested that 2D rolls undergo self-tuned nonlocal wavy instability ||13|| , which prevents 
their further growth. The nonlinear superposition of two sets of wavy rolls may result in 
the form of square patterns. However, this proposition is not analyzed so far. 

We present in this article a simple dynamical system, which describes nonlinear inter- 
action between mutually perpendicular sets of wavy rolls, and captures the mechanism of 
selection of square patterns in zero Prandtl number Boussinesq fluid. We show that the 
generation of vertical vorticity is important in addition to higher order modes to provide 
nonlinear coupling among mutually perpendicular sets of rolls. This is qualitatively different 
from the case of high Prandtl number convection |TJJ where nonlinear interaction of two sets 
of straight rolls may give rise to square patterns even in absence of vertical vorticity. The 
mutually perpendicular sets of wavy rolls interact through distortions in the vertical velocity 
as well as the vertical vorticity. The nonlinear interaction give rise to complex convective 



2 



patterns. The sequence of wavy stripe along x— axis, square patterns, and y— axis is chaotic. 
The generation of square patterns from wavy rolls is via forward Hopf bifurcation. 

We consider a thin horizontal layer of fluid of thickness d, confined between two con- 
ducting boundaries, and heated underneath. The fluid motion is assumed to be governed 
by zero Prandtl number Boussinesq equations [|I|, |T^], which may be put in the following 
dimensionless form: 

d t V 2 v 3 = W A v 3 + RW H 2 9 

- e 3 .V x [(w.V)v- (v.V)w] (1) 
d t uj 3 = V 2 u 3 + [(u.V)v 3 - (v.VV 3 ] (2) 
V 2 9 = -v 3 (3) 

where v = (vi,V2,v 3 ) is the velocity field, 9 the deviation from the conductive temperature 
profile, and uj = V x v = (ui, u>2, w 3 ) the vorticity field. The Rayleigh number R is defined as 
R = gi^pM ; where a is the coefficient of thermal expansion of the fluid, g the acceleration 
due to gravity. The unit vector e 3 is directed vertically upward, which is assumed to be 
the positive direction of X3-axis. The boundary conditions at the stress-free conducting flat 
surfaces imply 9 = v 3 = d 33 v 3 = at x 3 = 0, 1. and Vjj 2 = dn + 822 is the horizontal 
Laplacian. 

We employ the standard Galerkin procedure to describe the convection patterns in the 
form two sets mutually perpendicular sets of wavy rolls, and the patterns resulting due 
to their nonlinear superposition. The spatial dependence of the vertical velocity and the 
vertical vorticity are expanded in Fourier series, which is compatible with the stress-free flat 
conducting boundaries and periodic boundary conditions in the horizontal plane. As DNS 
shows either standing waves or stationary patterns, all time-dependent Fourier amplitudes 
are set to be real. The expansion for all the fields are truncated to describe convective 
structures in the form of straight cylindrical (2D) rolls, wavy (3D) rolls, and patterns arising 
due to their nonlinear superposition. The vertical velocity v 3 and the vertical vorticity uo 3 
then read as 

^3 = [Wioi(i) cos (kxi) + W m i(t) cos (kx 2 )] sin (nx 3 ) 
+ W / 2ii(^) cos (2kxi) cos (kxz) sin (tcx 3 ) 
+ Wi2i(t) cos (kxi) cos (2kx 2 ) sin (ttx 3 ) 
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+ Win (t) sin {kx\) sin (kx 2 ) sin (710:3) 

+ W\i2{t) cos (kxi) cos (/cx 2 ) sin (27rx 3 ) + .... 

u 3 = Z 010 (t) cos (kx 2 ) + Z 100 (t) cos (kxi) 

+ Zm(t) cos (fcxi) cos (fca^) cos (7nr 3 ) 

+ [Z 2 oi (0 cos (2kxi) + Z 02 i (t) cos (2kx 2 )] cos (7nr 3 ) 

+ [Zio2(0 cos (kx±) + Z m2 {t) cos (fcx 2 )] cos (2^x3) 

+ Z 2 io(i) cos (2kx\) cos (/cx 2 ) 

+ Zi2o(t) cos (fcxi) cos (2/1x2) + •••• 



(4) 



(5) 



The temperature field is slaved in the limit of zero Prandtl number, and it can be computed 
from Eq. 3. The horizontal components t>i,t> 2 of the velocity field and the same (001,002) 
of the vorticity field are then computed by their solenoidal character. We now project the 
hydrodynamical equations 1 — 3 into these modes to get the following dynamical system 

k = 3 Ar-l)X+ 1 - 
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and 77 = Z\\\. All roman characters stand for velocity modes, and greek letters stand for 
vorticity modes. 

There is no fixed points in this model as growing straight rolls are exact solutions of the 
problem. We now numerically integrate this dynamical system to investigate time dependent 
solutions. For each value of r, we start integration with randomly chosen initial conditions. 
We integrate for long periods (~ 100—200 times the viscous diffusive time scale) to ignore the 
transient effects. We increase r in small steps Ar (= 0.01) and repeat the above mentioned 
procedure. We get chaotic solutions for 1.11 < r < 1.42 for k x = k y = k c = which 



is in qualitative agreement with the DNS of Thual |12j| . For r < 1.11 we observe wavy rolls 
oscillating chaotically, and for r > 1.11 we observe chaotic sequence of wavy rolls, squares, 
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FIG. 1: Four different patterns varying chaotically in time for r = 1.11: (a) wavy rolls along x-axis, 
(b) square pattern, (c) patchwork quilt pattern , and (d) wavy rolls along y-axis. 



and asymmetric squares. For r > 1.42 and k x = k y = k c , the model is not good enough as 
the anharmonicity in various fields is very high in the limit of zero Prandtl number. This is 
natural as there is no additional heat flux across the fluid layer due to convective motion. 
The extra energy is consumed by the internal dynamics. 

Figure 1 shows competition of various patterns for r = 1.11. Two mutually perpendicular 
sets of wavy rolls compete with each other. This leads to square as well as patchwork quilt 
patterns. The sequence of occurrence of these patterns is not periodic but chaotic. Notice 
that the square pattern (see Fig. 1) consists of two sets of squares. A small (big) square 
has four big (small) squares as its nearest neighbors (Fig. 1 b). This feature of the square 
pattern is qualitatively new for convection in Boussinesq fluid. 

The top row of Fig. 2 (starting from left) shows time evolution of velocity modes W\oi, 
Won and W112 respectively after integration of the model for a period more than 200 x 
viscous time scale. The signal is chaotic. The second and third rows of Fig. 2 show various 
projections of the phase plots in fifteen dimensional phase space. Figure 3 shows time 
variation of the spatially average energy E av — ^ < v\ + v\ + v\ > Xl x 2 x 3 ■ The averaged 
energy shows relaxation oscillation in irregular way. This is also a common feature of 



DNS pj. 

We have presented in this paper a simple model which explains the mechanism responsible 
for generation of square patterns in zero Prandtl number limit in Boussinesq fluids. As the 
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FIG. 2: The top row (starting from left) shows chaotic evolution of Wioi, Won, and W\\2 for 
r = 1.11. The second row shows projections of phase space plot in Won — Wioij Wxi% — Wloi, 
and W\\2 — Won planes respectively. The third row contains projections of phase space plot in 
Zioo - Z ow , Z 010 - Wioi, and W n - Z wo planes respectively. 




FIG. 3: Time variation of spatially averaged energy in a unit cell. 



Rayleigh number is increased, one set of wavy rolls oscillating chaotically become unstable. 
Another set of wavy rolls are generated in a direction perpendicular to the former one. The 
nonlinear superposition of these wavy rolls give rise to squares and other complex patterns. 
The vertical vorticity modes with nonzero mean in vertical direction are responsible for the 
wavy nature of rolls at very close to the instability onset in case of stress-free bounding 
surfaces. These modes stop the unlimited growth of the 2D rolls. The DNS [O also showed 
wavy rolls rather than straight rolls. The nonlinear modes depending on both the horizontal 
coordinates facilitate the exchange of energy between two sets of wavy rolls. 
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